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Abstract 

We develop series representations for the Hurwitz and Riemann zeta func- 
tions in terms of generalized Bernoulli numbers (Norlund polynomials), that 
give the analytic continuation of these functions to the entire complex plane. 
Special cases yield series representations of a wide variety of special functions 
and numbers, including log Gamma, the digamma, and polygamma functions. 
A further byproduct is that C( n ) values emerge as nonlinear Euler sums in terms 
of generalized harmonic numbers. We additionally obtain series and integral 
representations of the first Stieltjes constant 71(0). The presentation unifies 
some earlier results. 
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Introduction and statement of results 



Let C( s ) a ) be the Hurwitz zeta function and £(s) = £(s, 1) the Riemann zeta 
function [TOJ [T3J [TS] . We develop series representations for these functions in terms of 
coefficients with the generalized Bernoulli numbers and these have several conse- 
quences. We illustrate that our results also extend to Dirichlet L functions. As Corol- 
laries, we obtain series representations of the log Gamma function, the polygamma 
functions, and special numbers including the ordinary and generalized harmonic num- 
bers. Specific application of our series is to the Stieltjes constants 7&(a), leading to 
new integral representations. We are able to unify some earlier results, specifically 
including the very recent ones of Rubinstein [13] . The Discussion section gives details 
as to how our framework subsumes that presentation. While an experimental math- 
ematics approach lead to a representation and demonstration of certain derivative 
values oc' k {l), this is a natural result within our approach. We then provide some 
other observations concerning the Norlund numbers implications of our series 

representation of ((s,a), and finish with selected concluding remarks. 

The generalized Bernoulli numbers are explicitly given by [33] 




(1.1) 



wherein a may be complex. On the right side, a factor with the Stirling number of 



the second kind S is evident, 
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B^ is a rational polynomial of degree n, with highest coefficient (—1/2)". The first 
few such are B^ a) = 1, B[ a) = -a/2, B { 2 a) = a(3a-l)/12, and B { 3 a) = -a 2 (a- l)/8. 
For these Norlund polynomials, a = 1 is a simple root if n > 1 is odd, and a = is 
a simple root if n > is even, and a double root if n > 1 is odd. 

The generalized Bernoulli polynomials (e.g., [ElEEE]) are given by 



^ ) (^) = EQ4 Q) ^- fc , (1-3) 
B}^'{x) = (— l) n B^(a — x), and they enter the generating function 

e - = $>(«)(*)-, |*|<2tt. (1.4) 

e i/ n=0 ft. 

As usual, B$p(x) = B n (x) and B^'(0) = B n (0) = B n denote the Bernoulli polyno- 
mials and numbers, respectively. 

Throughout we write a = Re s. In the following T is the Gamma function, 
(a)j = r(a + j)/T(a) is the Pochhammer symbol, ip = T'/T is the digamma function, 
7 = — ^(1) is the Euler constant, ip^ is the polygamma function, and p F g is the 
generalized hypergeometric function (e.g., [T| |2| |9| IT9]). 

The Hurwitz zeta function, defined by ((s,a) = J2^=o( n + a ) s for a > 1 and 
Re a > extends to a meromorphic function in the entire complex s-plane. In the 
Laurent expansion 

1 00 (—i)n 

C{s, a) = — - + £ { —T^(a)(s - IT, (1.5) 

S 1 n =0 n - 

7„(a) are the Stieltjes constants [H IU [T7J [20], and by convention one takes 7^ = 



We have 

Proposition 1. Let s ^ 1 and Re a > 0. Then we have 

U ' J r(s-l)^ fc! ( s + fc-l)r(s + fc + a-l)' 1 J 



In particular, 



11 i oo (_-t\k f>( 

C(s) = — - + - Viii f* fi7) 



Equivalent ly to (1.6) we have for Re a > 0, 



00 f— 

T(s)[s((s + l,a)-a~ s } = T(a)J2 { " ' 



fc-i 



fc! 

fc=0 



(s - l)Bi s+k - 1] sB { k s+k) 



r(s + fc) 

r(s + fc + a)' 
(1.8) 



s + fc — 1 s + fc 

Corollary 1. For integers n > 1 we have 

,<•>(*) = <_ ir ^)g<^|^±^. (, 9) 

This representation holds due to the relation ^ n \x) = (— l) n+1 n!C(^ + l,x). In 
particular, for the trigamma function we recover 

^' {x) = T(x) £ (fc + i)r(x + fc + i) = § (fc + i) (* + fc)R? (L10) 

and the functional equation tp'(x + 1) = ip'(x) — 1/x 2 : 

In the case of the tetragamma function, the functional equation follows from 

oo r,| TT 

rix + 1) = nx) - 2r(x) e T{x+k k +2y ( L12 ) 
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where Hk is the fcth harmonic number, and the latter sum is given by 
^ fc! * 1 = ~ 1 ~ fc! 



1 oo 

E 



1 



(1.13) 



xf^Tix + i + l) x 3 T(x)' 
In connection with (1.10)-(1.13), some elementary relations have been relegated to 

Appendix A. 

We recover a well known value in the next result. 

Corollary 2. We have 

<Z(0,a) = ±-a = -B 1 (a), (1.14) 

where Bi is the first Bernoulli polynomial. 

We also have series representations for the digamma and log Gamma functions, 
and for the first Stieltjes constant. 
Proposition 2. Let Re a > 0. Then (a) 



-l)n I 



*«) -■».= ,£ „ (a)n 



(b) for Re x > -1, 

Inr(x+l) = 5:m— B«(l) 



(1.15) 



fc=0 



fc! 



x r(fc - i)r(x + 1) 



fc- i fc 



r(fc + x) 



;i.i6) 



and (c 
1 



-1 



,n-l 



n— 1 



;i.i7) 



We emphasize that the sum on the right side of (1.16) is well defined. For fc near 1 we 
have T(fc — 1) = l/(fc — 1) — 7 + 0(fc — 1) and for fc near we have T(fc — 1) = — l/fc + 
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7 — 1 + 0(k). Then the quantity on the right side in brackets is —x[ip(x) + 7] + x — 1 
for k = and tp(x + 1) — x + 7 for k — 1. 
Proposition 3. (a) For Re a > we have 

7i(o) = \W{a) - ^ 2 (a)] - - f 1 v 4 F 3 (1, 1, l,v + 1; 2, 2, a + 1; l)cfo, (1.18) 
z a Jo 



(b) we have 



and 



and (c) 



7i = y + ^ (7^) + - ^'0*0]) <fo, (1-19) 



7i 



7T 2 1 



--/ 

2 Jo 



2- + ^'(a:) -^ 2 (x) 

3v 



dx, 



(1.20) 



72(a) 



-^ 3 (a) + 3tp(a)ip'(a) - ip"(a)] + 2T(a) J2 

00 ( — 



1 



+2r(a) 2 



fc= "! A; 3 T(A; + o) 
Corollary 3. We have 



W T(k + a) 

) 

'-l+-fk + k^(k)-ktp(k + a)]. 



ds 



J s=l 



(1.21) 



00 f(h\ 00 ( ■\\m m— 2 

2[C(2) - 7i] = 7 2 + 1 + 2 E(-l) fe ^T " E ^37 E C(*)C(m " (1-22) 



fc=2 



m=4 



fc=2 



A sum of (1.22) has many alternative forms, a few of which are collected in the 
following. We let T(x ) y) be the incomplete Gamma function. 
Corollary 4. 

gnr c(A;+1) = / b+m^+int]-^ 
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i=i 



J \ 3 



OO 



[7 + V(l + e -t )]dt= y [7 + V(e~*) +e*]dt (1.23) 

This sum, with approximate numerical value 1.25774688694, has been encountered 
before in analytic number theory [6] in the form Y^T=i + + !)]• The latter 

reference gives several other representations of this sum. 

A recursion for the derivatives of the Norlund polynomials is given in the following. 
Proposition 4. For n > 0, 

With a (null) convention for the sum, this includes the initial cases d a B ( a) = and 
d a B[ a) = -1/2. 

Proposition 1 can be applied to Dirichlet L-functions. As a first quick example, 
the Dirichlet L-function defined by 

oo (_\)n 

corresponding to quadratic characters modulo 4, can be expressed as 

L(s) = 4- s [((s, 1/4) - C(5, 3/4)] = 1 + 4" s [C(s, 5/4) - ((s, 3/4)]. (1.26) 

In particular, we have for nonnegative integers m the special values 

L(2m + 1) = - 2(2w ; + 1)! ^ + i(l/4). (1.27) 

Generally the values of L at odd or even integer argument may be expressed in terms 
of Euler or Bernoulli polynomials at rational argument and these in turn expressed in 
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terms of the Hurwitz zeta function. Therefore we may in this way obtain many other 
computable series representations for L(2m) and L(2m+1). These include the special 
cases of L(l) = tt/4, L(2) = G ~ 0.91596559, Catalan's constant, and L(3) = vr 3 /32. 

In general, Dirichlet L functions may be written as a combination of Hurwitz zeta 
functions. For instance, for x a principal (nonprincipal) character modulo m and 
a > 1 (a > 0) we have 




(1.28) 



Proof of Propositions 



Proposition 1. We will use the generating function [12] (p. 147) 




x\ < 1, 



(2.1) 



valid for complex z. We have for o > 1 





1 



E 



oo 



(-i) fe 



B(s + k — l,a) 



(2.2) 



fc=0 



fc! (s + k-1) 



where the Beta function factor 



B(s + k — l,a) 



T(a)T(s + k - 1) 
r(s + k + a- 1) ' 



(2.3) 
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With B^ +k ^ being polynomials in s, (2.2) gives the analytic continuation of ((s, a) 
to the whole complex plane, completing the Proposition. 
(1.8) is based upon the use of (2.1) and the integral 



/' 

Jo 



1 1 
ln(l - y) y 



(1 - yf- 1 In* f - J- ] dy = r(s)[s((s + 1, a) - a~ s }. (2.4) 



Based upon (2.11) below for a~T(s), we see the equivalence of (1.7) and (1.5). For 
(1.7) we may note the particular value [B^" 1 ^ jik — 1)] fc=1 = — J^tdt = —1/2. 

Corollary 2. We recall that B^ a) = 1 and b[ s) /s = -1/2. When s ->• in (1.6), 
only the k = and k = 1 terms contribute, as l/T(s — 1) = — s + (1 — 7)s 2 + 0(s 3 ). 
We obtain C(0, a) = 1 - a - 1/2 = 1/2 - a. 

Remarks. With B^(x) the generalized Bernoulli polynomial for v < n, 

j/l d n ~ v 7/' d n ~ v 
Bi n+1 \x) = -- (x - l)(x - 2) • • • (x - n) = (-1)"-^ (1 - x) n . (2.5) 

We have = 5^(0). It follows as a very special case of (2.5) that B^ 1 ^ = 
(—l) h k\. This latter relation allows the recovery of the series £(2) = Y^=i l/^ 2 from 
(1.7). 

As also follows from (2.5), B ( k k+2) = (-l) k k\H k+1 , where H n = £Li 1/fc is the 
nth harmonic number. With these values, we have from (1.7) an Euler series (e.g., 
[8J, Appendix B) for C(3): 

-1 OO TT OO TT 

Integral arguments of the digamma and polygamma functions are directly related 
to the harmonic H n and generalized harmonic Hj?' numbers: H n = tp(n+l) —-0(1) = 
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i>(n + 1) + 7, 



H n ] = S [^Kn + 1) - ^(l)] > ( 2 - 7 ) 



Therefore Corollary 1 gives another form of these special numbers. The generalized 
harmonic numbers are given by 



n -1 

j= iJ 



As also follows from (2.5), we have the relation 

B^) = (-iyn\[Hl +2 -H%]. (2.9) 
Proceeding as for the Proposition, we have for Re a > 0, 

roo pi 

a~ s T(s)= x s ~ 1 e~ ax dx= [-\n{l-t)] s - l {l-t) a - l dt, (2.10) 
Jo Jo 

so that using (2.1) we have 

K J v ; ^ k\ (s + k - 1) T(s + k + a) v ; 

We have found that (1.9) and (1.15) are given in Norlund's book [l2j in a section 
dealing with numerical differentiation and integration (pp. 243 and 244, respectively). 
Proposition 2. (a) We introduce the constants ([5], Proposition 11) 

Pn+i = r / {-x) n dx = - 2^ , 2.12 

n! Jo n\ i~ x k + 1 

where s(k,£) is the Stirling number of the first kind. These constants enter the 

generating function 

oo 1 I 

Y,Pn + iz n ~ l = - + — T, |*| <1. (2.13) 
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Lemma 1. We have 



(2.14) 



We have 

[\x) k dx= (\l-y) k dy={-l) k B < £ ) = [\k - y)(l - y) k . x dy 
Jo Jo Jo 

= - 1 2/(1 - y)h-xdy + k [ (1 - y) k -idy, (2.15) 
Jo Jo 

giving (— l) k B^ = —k\pk + i + k(—l) k ~ 1 B < £i l 1 \ from which the Lemma follows. 

Part (a) of the Proposition then immediately follows from Proposition 5(a) of [7J. 
Here we give direct verifications: 



E 



-i) n i 

i n (a) n 



l l r 1 



n (a n Jo 



n=l 



[(x) n - n(x) n - 1 }dx 



°° 1 1 f 1 f 1 

= ——— / (x — l) n dx = / [if) (a) — if) (a — x + l)\dx = if) (a) — In a. 
^ n [a) n Jo Jo 

Otherwise, we may employ 
1 r(a) 



(2.16) 



{a) n T(a + n) T(n) 



1 r°° 

— / e- at (l - e-y^dt, Re a > 0, 

,71) Jo 



(2.17) 



to write 



E 



:-i) n i 

i n (a) n 



4 n) + nSfc^ 
1 



1 /"OO (p* p xt\ 

at{ 'dtdx 

o Jo 



-at 



e"* - 1 t 



(1-e* 
(it = VK ) — hi a- 



(2.18) 



The latter integral representation is a standard one (e.g., [9], p. 943). 
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(b) From (2.1) and B ( ^ +1 \l) = (l - J) it follows that 



C \ z oo k 

w+v) =£^* +1,(1) ' (2 ' 19) 

We use the integral representation [TT] (p. 343) 

(2.20) 

with the last integral given by the Beta function £?(&; — 1, x+ 1), and part (b) follows. 

For (c) we may use Lemma 1 and Proposition 5(b) of [7], or else apply the integral 
representation (2.4). 

Remarks. It is easy to show that 

= E(-l) fc TZT' (2 ' 21) 
where s(n, n) = 1 and s(n, 0) = 5 n Q in terms of the Kronecker symbol 5jk- 

From m (p- 147), (k - l)flf(l) = -S?" 15 = (A; - ^(-l)^ 1 Jo *(1 " t)t-idt, 
so that for the coefficients of part (b) we have B^ k \l) = (— l)*" 1 Jq 1 t(l — t)k-idt = 

In fact we have P2] (p. 148) 



n 

-l) n (l-^)n = E(l> r ^ 1) - ( 2 - 22 ) 

r=0 



r. 



Therefore, performing the integrations, we obtain 

,(n+l) 



/ (x) n dx = (-irE U^-, (2.23) 

^ w (r + 1) 
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for n > 0, 

(-!)• / o ,(1 - = = (-1)" g (.,,) ^ (2-24) 

and 



( -ir/ o H^ = (-irEjj(2 r+1 -i)^. (2.25) 

We record the latter result in the following form for the constants of (2.13). 
Lemma 2. 

P.. +1 -^E(:j(2" + '-l)^ y . (2.26) 

From differentiating the generating function (1.4) there is the relation B^ +1 ^ = 
^B^-kB^. At a = A; -1 this yields k\p k+l = (-l)*" 1 ^* _1) /(A;-1). According 
to the relations above (2.22), this expression is equivalent to (2.12). 

Comparing (2.1) with the generating function 

oo k 

M* + l)] B = n!E*(M)fcP ( 2 - 27 ) 

we have 

(m-lV B (m) 

^' t-m^v (2 - 28) 

A generalization of Stirling numbers of the first kind in the first argument is presented 
in Appendix B. 

Proposition 3. With B^ = 1, Proposition 1 gives 

r(q) i r(q) ^ (-i) fc iff 4 *- 1 * r( 8 + fc-i) 

U ' J r(s + a- 1) (s- 1) r(s-l)^ k\ (s + k — 1) T(s + k + a - 1)' 

(2.29) 
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We then expand both sides about s = 1, using for instance 



V{ s , } (a - 1) + 7 (« - I) 2 + [y - 12 ) {S - 1)3 + ° [{S ~ 1)4] ' (2 - 30G) 



and 



r(s ^ ( " } _ T) = 1 - «<<)(» " 1) + \W» ~ ^(»)](» " I) 2 



+ 6 



-^ 3 (a) + 3^(a)V'(a) - V"(a) (s - I) 3 + 0[{s - l) 4 ]. 



(2.306) 



From the 0(s — 1) terms we recover the known relation 70(a) = —ip(a). In addition 
we have 



C(s,a) + 



{8-1)* 

Now [12J (p. 147) 



S->1 + 



1 00 ( — ~\) k r(*) 

-7i(«) = 2^ 2 («)-^'H] + r(a)E "' J/ 



, = ; A; 2 r(* + a)" 
(2.31) 



flf = (-l) fc [\l-t) k dt = (-l) k [\t) k dt. 
Jo Jo 



(2.32) 



1 (t) k 



1 



-t 4 F 3 (l,l,l,t + l;2,2,a + l;l), (2.33) 



We easily have 

^fc 2 r(A; + a) T(a + 1)' 

and part (a) readily follows. 

Part (b) then follows at a = 1. For (1.20), in light of ip(x) = — 1/x — 7 + (7r 2 /6)x + 
0(x 2 ) as x — > 0, we write the term ^(x) = ^[^(x) + 1/x — 1/x]. We then employ 
the integral 



L 



ij){x) + - 
x. 



dx= [ il>{x + l)dx = lnT(x + 1)|J = 0. (2.34) 
Jo 
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Part (c) proceeds similarly to (a), also using 



D h — Du 



dB 



(s+k-1) 



ds 



{*-!) + 



d 2 B { ° +k - l) 
ds 2 



Expression (1.21) follows from 



CM- 



72(a)- 



{s-lf + 0[{s-lf\. 

(2.35) 



(2.36) 



s^l+ 



{s-lf_ 

Remark. We have the Laurent expansion about s — 1 of the logarithmic derivative 
of the zeta function, 



(2.37) 



where rj = —7. From Proposition 2 or Corollary 3 we then obtain corollary expres- 
sions for rji = 7 2 + 27! . 

Corollary 3. We use Proposition 3(b) along with the expansions (e.g., [9], p. 944 
or p, p. 259) 



V(x + 1) = ${x) + - = -7 + E(-!) fc C(fc) 



.r 



fc=2 



1 00 

^) = - + E(- 1 ) fe+lfc c(^ + i) 
x fc=i 



;r 



jfc-l 



|x| < 1, (2.38) 
(2.39) 



and 



1 00 00 

^{x) = - + 2 1 - + 7 2 - 2J2(-I) k ((k)x k ~ 2 - 2 7 E(-l) fc C(^ 

X X fc=2 fc=2 



00 m— 2 

+ E £(-i) m C(£0C(m-£0z 

m=4 fc=2 



ra-2 



(2.40) 
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We insert these into 

2[C(2) - 71] = f [2^ + tf(x) - ^ 2 {x)\ dx, (2.41) 

J . X 

integrate term by term, Abel sum the series 

00 

lim ]T(-l) fc+1 C(A; + l)x k = 1, (2.42) 
k=i 

and use the easily proved sum I]£L 2 ( — l) fc CWA = !■> giving the Corollary. 

Corollary 4- The first expression follows by using a standard integral represen- 
tation for the zeta function. The next follows by use of a geometric series, and 
the integral representation with the digamma function follows by using the Laplace 
transform representation of 1/k. 

Proposition 4- From the generating function (1.4) we have 

e ln h^r =E— Ir^rr. 1*1 < 27r - ( 2 - 43 ) 



e z — 1/ \e z — 1/ ~[ da n\ 



In order to expand the log factor, we first note that 

/ z \ 1 1 1 1 ~ „ z n . 

^— =--1-^— = --l--]T£ n -, 2.44 
\e 2 — I J z e z — 1 z z n^o n - 



dz 
yielding 



d f z \ , ^ 5,,+! 
m 



Upon integrating, 



I ( Z \ \^ -5n+l^" +1 _^ _ -5n+l^ W+1 

n ^*-iJ~~ * ^(n + l)(n + l)! " 2 ^ (n + l)(ra + 1)!' " 
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Then at x = in (2.43) we have 





z 


d 1 




00 


00 



da rcJ ^ n 

n=l ^ "'■ n=0 

1 oo n+ 

= __V J B (a) - V V ±Jn '^ m+l ^ (2 471 

2£S » n! fe,fen!(m+l)!(m+l)- 1 ' J 

Reordering the double sum, 

^ 9a n! 2^ » n\ ^^ n !(* - n + 1)! (/- n + 1)' l ' j 
from which the Proposition follows. 

Remarks. It is seen that the recursion (1.24) also applies for x 7^ 0. Since Bm-\ = 
for n > 1, the sum of (2.44) may be written as \ Em=i m(2m)T • From (2.43) we have 
similarly for higher order derivatives 

: • / z \ 00 d j B (a Ux) z n 



implying 



u^m (-0) 



, 9a- 7 / n n\ 



Discussion 

Rubinstein [T3] developed several expansions for the Riemann zeta function using 
certain polynomials afc(s). That work employed the generating function 



ln(l-t)V _1 



V 1 J 

Comparing with (2.1) we identify 



5> fc (s)t fe , |t|<l. (3.1) 

k=0 



B (s+k-l) 



ai ( s) = (s _ 1 )^_^_. (3.2) 
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We have 

a' (s) - ^ + t¥ d ( 3 3) 

a *W- (]fc-l)!( a + Jfc-l)2 + fc! ( S + Jfc-l)da fe ' 1 ' 

where, by (1.1), and using the functional equation of the digamma function, 

d a Bj?> = [V>(n + a + 1) - ip(a)]B^ 

-£(:-:)(°^ 4 "^^£^C)^ (3 - 4) 

and ^(n + a + 1) — ^(a) = S"=o V( a + r )- We have precisely by (2.32) 

^"TET-r-lSijiw**' (3 - 5) 

that was found by indirect means in [14] . 
We also have for k > and m > 0, 

(-D fc 

= Sr^' 3.6a 

a 2m+2 (-2m - 1) = ^ 2m 2 " + ^, , (3-66) 

and 

« 2m+2 (-2m) = -(2m + l) ,.^' ■ (3.6c) 

(2m + 2)! 

By using the generating function (1.5), one easily sees that B^ 1 ^ = l/(n + 1), so that 
from (3.2), a fc (-Jfe) = (A; + giving (3.6a). 

On the Norlund numbers 

These numbers have the known asymptotic form Bjfi ~ (— \) n n\/ Inn as n — > oo. 
From (2.32) it is easy to see how a standard generating function for them arises: 

oo r>(n) r i 

y -n- z n = / (i + z )-*dt = -_ (4.i) 

to n\ Jo 1 ] (l + z)ln(l + *) 1 ; 



From (2.32) and the relation above (2.22) for B^ T 1 - ) it is also clear how to develop 
expressions given in [12J (p. 244): 

h ^ = ££W (4 . 2) 



In 



x J ^ (x + 1) 
and 

's + lW 1 1 (-l) r ^ r " 1) 1 f43l 

x J x 2x(x + l) x^ 2 (r-1) (x + l) r ' l ' J 

We note that differentiation with respect to x of these formulas gives further sum 
identities for generalized Bernoulli numbers. 

Given the generating function (2.19) and that for harmonic numbers, 

H k+ iz k = f- f , 4.4 

one suspects various relations between generalized Bernoulli numbers and harmonic 
numbers. As an example, we present the following. Although this relation could be 
proven with generating functions, we give a proof employing special function theory. 
Proposition 5. For integers n > 0, 

n I -r \n—r 

Er^tt + i = n + l. (4.5) 
r=0 (n-r)! 

Proof. We begin by reordering sums and applying Chu-Vandermonde summation, 

n I i \n— r n ( i\n-r r+1 -r n+1 i n / i \n—r 

\ " I 1 J r>(n-r) ff I r>(n-r) 1 _ 1 \ " I p}(n-r) 

n+1 in i .1 



^1 r 1 r(t + n-l + 2) 

£t r(t + i)r(n-i + 2) 
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= -!—[ 1 (lL— )dt = n + l. (4.6) 

n+ 1 Jo \dtB{n + l,t+ 1) J K J 

Above, we manipulated a terminating 3 F 2 function at unit argument, since 

^i r(t + ra-l + 2) " (l) 2 e T(t + n-£ + 1) z e+1 
T(n-£ + 2) Z ~f^(2) e T(n-£ + l) t\ 

^(2), T(n+1) (-t-n)* T(n + 1) 21 ; 

(4.7) 



Because 



we have 



3 F 2 (1, 1, a; 2, b; 1) = {^zjj [^(6 - 1) - ^(6 - a)], (4.8) 



r(* + n + l) n/1 , rt r(t + n + 1) (t + n + 1) . ,. „ N , . 

3 F 2 (l, I, -n; 2, I) = r ( n + l ) („ +l ) l^(~ t ~"~ l )~^'(~ t )l 

^ r( r(t'+2) 2> Wt + " + 2> "^ i + l)l ' (4 ' 9) 

by recalling that ^(l — x) — ip(x) = ncot^nx). 

Other relations from (1.6) 

From the representation (1.6) we have 

d a ((s, a) = -s((s + 1, a) = ip(a)((s, a) 

T(a) ™(-l) k B[ s+k - 1] V(s + k-l) , 

E J 77h— TTTv^ , , „_V (* + * + a - 1). (5.1) 



r(s-l)^ fc! (s + k — l) r(s + A; + a — 1) 
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We also have the special case 

H S '2j-r( S -l)^ k\ (s + k-l)T(s + k-l/2)- {2 1)C(S) - (5 ' 2) 

It would be desirable to otherwise have proofs of the right-most equalities in (4.1) 
and (5.2). The duplication formula for the Gamma function may be useful in this 
regard for (5.2). 

The Bernoulli polynomials satisfy the multiplication formula 

m-l / k \ 
B n {mx) = ffl"" 1 Yl B n [x + - , (5.3) 

that is easily verified with the generating function (1.4) with a — 1. One may ask 
whether there is a generalization to a formula such as 



B^(mx)=m n - a f 1 (n,a)Y: 



m—1 

(5.4) 



( x + ^.)+f 2 ( x , n ,a) 
\ m 



k=0 

where fi and fi are such that — \ and f 2 — for a — 1. 

Summary remarks 

The series (1.6) and (1.7) converge everywhere in the complex plane and include as 
special cases the relations £(1 — m,a) = —B m {a)/m and the trivial zeros £(— 2n) = 0. 
The treatment has been sufficiently general to subsume expansion of the zeta function 
at integer argument in terms of Stirling numbers. We may emphasize, among other 
features, that nonlinear Euler sums for ((n) naturally emerge as special cases. As an 
example, we have from (2.9) 
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in terms of generalized harmonic numbers H%' . 

Proposition 3 may be the only known integral representation for the first Stieltjes 
constant 7x in terms of the digamma and trigamma functions. 

It seems that the second line of (3.4) for the derivatives d a B^' could have an 
alternative, more compact form, possibly involving B^ and B n values. Such an 
expression would be very convenient in further developments for the Stieltjes and 
other constants. For instance, it could be immediately applied to Proposition 3(c). 

As mentioned in connection with (5.4), further investigation of the Norlund poly- 
nomials appears to be in order. 

Appendix A: Beta function-based relations 

We have for Re x > and n > an integer, the special case of the Beta function 

rl 77 1 77 ! 

B(x,n+1)= t*-\l-t) n dt= =—— , (Al) 

J o x(x + I) ■ ■ ■ (x + n) [x) n+ i 

that may be verified by induction. So as relates to (1.10)-(1.12), we have 

V _^ i = V - / t x ~ l (l - t) k dt 

t (Hl)(x) w ^ o (k + l)J 1 ; 

= f 1 l?—]ntdt = ip'(x). (A.2) 
Jo t — 1 

Here, we have recalled (e.g. by [9], p. 943) 

ij)V)( z \ = f 1 1 l n i t dt, Re z > 0. (A3) 

Jo t — 1 

In regard to (1.12), we recall the generating function with harmonic numbers 

£^=5^ = ^, (A4) 
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where LL, is the polylogarithm function, and Hq = 0, and its integrated form, 

OO TT 1 

E^ fe+1 = ^n 2 (l-,). (A5) 

Then from (A.l) we have 

OO 1,177 OO TT r \ 

1 /-l 1 poo 1 

= - / t x - 1 \n 2 tdt=- / uV^cfo = — . (A6) 

2 Jo 2 Jo x 3 K J 

By binomially expanding the integrand of (A.l) and/or by differentiating with 
respect to x we obtain other relations. For instance, we have 



d x B(x,n+l)= / **- 1 (l-t)»lntdt= rT — [^-^+1)]=-^-^ 

(A.7) 



Alternative points of view of B(x,n + 1) = Ylt=o [ x +e) (") are m t erms °f partial 
fractions or as divided differences of 1/x. 

Appendix B: Generalized Stirling numbers of the first kind 

We let s(j, k) denote the Stirling numbers of the first kind with integer arguments. 
We proceed to generalize these with the first argument complex. In the following if) 
again denotes the digamma function, the polygamma functions, and = 
T(z + k)/T(z) the Pochhammer symbol. 

It is shown in [5] (Lemma 1) that 



ds 



'-i) j+i £\s(j + m + i). (B.r 
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From this we make the extension to A G C 



k-l 



(B.2) 



(k - 1)! \d S/ 

We then obviously have agreement with the usual Stirling numbers of the first kind 
when A is a nonnegative integer. When A is a positive integer, (s)a-i is a polynomial 
of degree A — 1 in s. Hence by the definition (B.2), s(A, k) — when k > A. Moreover, 
we show that these extended Stirling numbers satisfy the same recursion relation as 
their classical counterparts. 
Lemma Bl. We have 



s(A, k) = s(X - 1, k - 1) - (A - l)s(A - 1, k). 



(B.3) 



Proof. We first note (s)a-i = (s + A — 2)(s) A _ 2 so that from (B.2) 



s(X,k) = 



(-l) A+fc d_ 
(k - 1)! yds, 



fc-2 



:s) A -2 + (s + A-2) — (s) A _ 2 



s=l 



-1) 



A+fc 



(fc-1)! 



fc-2 



-l) A+fc (A; - 2)!a(A - l,k - 1) + ( £\ (s + A - 2)^(a) A _ 2 



s=l 

(5.4) 



By using the product rule we find for the last term 



fc-2 



|)~-' (S + A _ 2) ^ (8)A . 2 



'A;-2 N 



D*- 2 (s) A _ 2 | s=1 + (A - l)£>;- 1 (a) A -2|^i 



s=l 



= (-l) A+fc (A;-2)!^_ 3 J S (A-l,A;-l) + (-l) A+fc - 1 (A-l)(A;-l)! S (A-l,A;). (S.5) 
Then (B.4) becomes 

S ( A > fc ) = jr^K-l^ik-m^y.siX-l, A;-l) + (-l) A+fe - 1 (A-l)(A;-l)! S (A-l, fc)], 
(A; - 1)! 

(5.6) 
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and the Lemma follows. 
Remarks. We have 



ds 



(s) A = (s) A [^(s + A)-^(s)], 



(B.7) 



so that the 'higher' s(A, k) values may be obtained via Bell polynomials. We recall 
the harmonic H n and generalized harmonic numbers: H n — ip{n + 1) — = 
i>(n + 1) + 7, 



H (r) 



("1 



\r-l 



(r-1)! 

Our definition of s(A, /c) then extends as desired. For instance, we have 



(5.8) 



s(A,2) = (-1) A (A- 1)!^, s(A,3) = (-l)^ 1 ^^^ -lf&], (5.9a) 



and 



s(A, 3) = (-!)■ 



6 



[i/ A 3 _ 1 -3/J A _ 1 ^ 2) 1 + 2iJ^i]- 



r(3) 



(5.96) 



Here, it is understood that (A — 1)! = T(A). 
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